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ABSTRACT 

'66 : 

ON ' The process of photon splitting, becomes allowed in the presence of strong magnetic 

field. We calculate the influence of magnetized plasma on the photon splitting 
absorption coefficient. We calculate the refraction coefficients and polarization vectors 
^ | ■ with the inclusion of the vacuum terms at an arbitrary value of the magnetic field, and 

, then find photon splitting matrix element taking into account the terms that vanish in 

vacuum, but may be nonzero in the presence of plasma. We find the photon splitting 
| rate in plasma with the density typical for a neutron star star atmosphere and identify 

a region of the parameter space where the plasma effects are important. 

in : 

■ 1. Introduction 

oo : 

The possibility of photon splitting process in strong magnetic field around 
B c = 4.414 x 10 13 Gauss has been first noted in the early seventies (|Adler 1971 , Bialynicka-Birula 
fc Bialynicki-Birula 1970| ). Adler (1971) found that only one type of splitting is allowed, namely 

I — [ I + ||, where a perpendicularly polarized photon splits into two parallel polarized photons. 

This calculation took into account the vacuum dispersion and also the influence of matter 
with density that of a pulsar magnetosphere. The normal modes are polarized linearly in such 
• i-H ! conditions. The photon splitting absorption coefficient for this allowed reaction is 

b / B A 6 / huj \ 5 

.«?: * = » 12 bi sinS ) U?) ' (1) 

Since then a number of authors have derived the matrix elements, and absorption coefficients for 
this process using different approaches (Papanyan & Ritus 1972, Stoneham 1979) and also recently 
(Baier, Milshtein & Shaisultanov 1996, Baring & Harding 1997). These papers have settled the 
controversy sparked by Mentzel, Berg & Wunnner (1994) who suggested that the photon splitting 
rate may actually be a few orders of magnitude higher than previously thought. 

This exotic process has found a number of astrophysical applications. A natural environment 
where it may play a role is provided by neutron stars, since a number of them have magnetic 
fields in excess of 10 12 Gauss, and about a dozen of radio pulsars has spin down fields larger than 
10 13 Gauss, it has been found that photon splitting plays an important role in the formation 
of the gamma ray spectrum of PSR 1509-58, where it inhibits emission above 1 MeV ( |Harding" 



Baring fc Gonthier 19971) . Soft gamma-ray repeaters (SGR) form another class of objects where 
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photon splitting has been suggested to play an important role. There were three firmly established 
SGR's with and the fourth source was recently discovered ( [Kouveliotou 1997 ). All of them are 



characterized by short durations, below a second, and spectra with cut-offs around 30 keV. A 
number of arguments have been presented ( Thompson fe Duncan 1995| ) pointing that these sources 



are actually magnetars, i.e. neutron stars with magnetic fields reaching 10 15 Gauss, for which the 
main pool of energy is the magnetic field. Baring (1995) suggested that photon splitting may 
be responsible for the spectral cutoffs in the spectra of SGRs, invoking a possibility of photon 
splitting cascades. Thompson and Duncan (1995) showed that even in the absence of cascades 
photon splitting influences the shape of the spectra of SGRs. 

In this work we investigate the process of photon splitting not only in vacuum but also in 
the presence of matter with the density like that in a neutron star atmosphere, when the plasma 
dispersion may be important and for some range of propagation directions the normal modes are 
polarized circularly. The relevance of plasma effects can have a very substantial influence on the 
absorption coefficients. For example, Bulik and Miller (1997) have shown that plasma effects leads 
to a broad absorption like feature below ~ 10 keV for a large range of SGR emission models. In 
section 2 we analyze the polarization of the normal modes in the strongly magnetized plasma, in 
section 3 we calculate the photon splitting matrix elements and absorption coefficients, in section 
4 we present the results and we discuss them in section 5. 



2. Polarization of the normal modes, refraction coefficient 



We use the normal mode formalism to describe opacities in strong magnetic field. The 
dispersion equation is 



k x 



-'(J 



x E 



+ 



eE = 0, 



(2) 



where e is the electric permitivity tensor, fi is the magnetic permeability tensor and E is the electric 



Ruchadze 1975, and Meszaros 



field of the wave. Equation (|2|) has been discussed by e.g. |Ginzburg 

(1992). In general equation @ has three solutions: one representing the longitudinal plasma 
oscillations that do not propagate (Meszaros 1992, p. 69), and two perpendicular representing the 
electromagnetic waves. In the presence of the magnetic field there exist two nondegenerate wave 
solutions of equation (||), and therefore the medium is birefringent. We follow the convenient 
method introduced by Gnedin & Pavlov (1974) for the case of tenuous plasma to solve equation (|2|). 

The roots of the dispersion equation are 



n j 



ni ± 



(3) 



where j = 1,2 indicates which normal mode is selected. The refraction coefficients are given by 
the real part of rij, and the total absorption coefficients are proportional to the imaginary part of 
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rij; £j = (2lo I 'c)Im(nj) , and 

ni = 1 + \{e xx + e ra cos 2 9 + e Z2 sin 2 9 - e xz sin 26* 

(4) 

-Hxx - n yy cos 2 - nj sin 2 61 + nJ sin 20) , 
riL = \{txx ~ e yy cos 2 9 - e zz sin 2 9 + e xz sin 26* 

(5) 

+Vxx - V y y cos 2 9 - fi zz sin 2 6> + /j,J sin 20) , 
i 

n c = -ji^xy cos 9 + e xz sin 9). (6) 

Here /z" 1 is the ab component of the inverse fi tensor, and we use the system of coordinates with 
the magnetic field along the z-axis, and the wave vector in the yz plane at an angle 9 to the 
magnetic field. 

We describe the polarization of the normal modes by the position angle Xj between the 
major axis of the polarization ellipse and the projection of the magnetic field B on the plane 
perpendicular to the wave vector k, and the ellipticity V, whose modulus is equal to the ratio 
of the minor axis to the major axis of the polarization ellipse, and whose sign determines the 
direction of rotation of the electric field: 



v i = ' Tj exp ( 2 *^') = h • 

The polarization vectors can be described by a complex variable b, or by two real parameters q 
and p (Pavlov, Shibanov, & Yakovlev 1980) 

b = q + ip = — . (7) 

nc 

Let us consider a wave travelling in the direction described by 9 and 4> m the coordinate 
system with the magnetic field along the z-axis, and the x-axis chosen arbitrarily. The wave vector 
is k = k (cos 9, sin 9 sin (p, sin cos 0). In the rotating coordinate system e± = 2~ 1 / 2 (e x ± ie y ), 
eo = e z , the polarization vectors are 

4 = P+ 1 ^=C 7 -e^(^cos0± 1) 

(8) 

e o = CjKj sin 9 



where 

Kj=b\i + (-iy(i + b-i) 1/2 

and Cj = (1 + \Kj\ 2 )~ 1 / 2 . In general Kj are complex. 



(9) 
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2.1. Vacuum polarization effects 

We use the convention where we label the polarization state by the direction of the electric 
field vector of the wave in relation to the external magnetic field. Adler (1971) used a different 
convention using the magnetic, not electric field of the wave to define polarization. 

The refraction indices in strong magnetic field have been derived by Adler (1971). The result 



is 

-,vac 



nil i — -l 



1 - -sin 2 9A^{ujsin6,B) , (10) 



where the functions A^' 1 - are rather lengthy double integrals given by equation (51) in Adler 
(1971): 

A\\> ± (u,B) = = -^ ^exp(- S )/ exp^ilCMF^M (H) 
2ir Jo s z Jo 

and v = 2ts~ 1 — 1. The remaining functions j" ,_L are 



Bscosh(Bsv) Bsv sinh(Bsv ) coth(i?s) 2Bs[cosh(Bsv) — coshBs)] 
sinh(l?s) ' sinh(i?s) sinh 3 (i?s) 



ii Bscosh(Bsv) „ , . 

j" = -X . ; - Bscoth(Bs) 

sinh(Bs) 



9 sinh(i?sw) 
1 - v 2 +v- 



sinh{Bs) 



The function R is given by 



R{s,t) = \ 



t\ cosh(Bsv) — cosh(-Bs) 
2t { 1 I ~r 



s ) B sinh(i?s) 



(13) 



(14) 



The refraction indices ri^ = 1 + rj sin 2 9 have been calculated for hu; <C m e c 2 and arbitrary 
magnetic field by Tsai &; Erber (1974) In the low field limit, B <C B c , 

/ -,\ 14^ 2 13 . 

m(h) sa — 7r /i 4 

/||V 7 45 315 

... 8 o 379 , 4 

t7_l(/i) sa — h hr , 

' K J 45 5040 ' 

and when B B c 

nCO ~ ^-fr'inpfc), 

where /i = B/B c , and Li = 0.249.., and 7 = 0.577.. is the Euler's constant (Tsai & Erber 1974). 
Thus in the limit of strong magnetic field refraction for the parallel mode grows linearly with the 
field while that of the perpendicular mode is nearly constant. We present the dependence of the 
vacuum refraction coefficients as a function of the mangetic field in Figure |]. 

We have evaluated numerically the integrals of equations ( |l0| ) . Figure |2| shows the dependence 
of the refraction coefficients on the photon energy for a few magnetic fields. 
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2.2. Refraction coefficient in the system of plasma and vacuum 

Plasma in strong magnetic field can be described by the dielectric tensor 

e ab — €-ab + e a b ~ °ab i 

and the vacuum permeability tensor. The plasma dielectric tensor can be expressed as 



LO, 



tab = 5ab- 4 Hoi, where lo p 



/ 4vriVe 2 



1/2 



111 



is the plasma frequency and H a b is the plasma 



polarization tensor. The plasma polarization tensor is diagonal in the rotating coordinates, and 
for cold electron plasma is given by 



LO 



LO + auJB — ijr (*H + CtUJB 



(15) 



where a = —1, 0, +1, and j r = (2/3)(e 2 /mc 3 )w 2 is the radiative width, and we denote oot = uo — i^ T . 
Inserting equation (|D|) into equations (||) and (|6|), we obtain 



4 lo 



riL 



;i + cos 2 ^)- 



L0 t 



LOT — LO 



+ sin 2 6- 



8111 



LO, 



LO 



L0 t 



+ - (v4.ll (u sin 9, B) + A L (u sin 9,B)), (16) 



LOLOf LO 2 — L0 2 B 



L + L4 11 (w sin 0, j5) - A l (lo sin 0, 5)] 



and 



n c 



L0 B 



2 LO \Loi 



cos 6 . 



LO 



(17) 



(18) 



B , 



Inserting equations ([H]) and into equation (|7]), and neglecting terms proportional to j 2 we 
obtain 



sin B 
2 cos w 



1 - {AW{Losm9,B)- A^(Losm6,B))^ [ ^ - 1 



\ a; 



and 



P 



sin 



2 cos i 



LO 2 



1 + 2(All (w sin 6, B) - A 1 - (lo sin 6, B) ^ 



(19) 



(20) 



polarization vectors and refraction coefficients are determined by equations (|l6|) , (|l7|),and (|ilj|) 
combined with (||) and (g). 

The presence of matter influences refraction also because of the electron cyclotron resonance at 
LOb. These effects increase for propagation direction along the magnetic field, and are proportional 
to the matter density. For a detailed discussion see e.g. Meszaros (1992). 
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3. Photon splitting rate 

We will use two coordinate systems U, and U'. The z-axis in the system U lies along the 
magnetic field, and the wave vector of the initial photon is in the zx-plane. The z-axis in the 
system U' lies along k, the wave vector of the initial photon, and the magnetic field lies in the 
xz-plane. The system U is convenient for calculation of the matrix element M and the refraction 
indices rij, while U' will be used for integration. In U the wave vector of the initial photon is 

= feo(cos 0, sin ^o)- The vector k\ in U' is k\ = (cos 0', sin 9' sin (/)' , sin 9' cos (j)'), and in U it 
is k\ = k\ (cos 6*o cos 9' — sin 9$ sin 9' cos 4>' , sin 9' sin <j>' , sin 9$ cos 9' + cos 9q sin 9' cos 4> ) . Using the 
momentum conservation we obtain in U' k<i = k — k±. Polarization vector of a photon is given 
by e(k) = e(u,9,(p). We calculate the polarization vectors for each photon separately, given its 
frequency and the direction of propagation. 

The photon splitting absorption coefficient can be found by integrating the S-matrix element 
over the phase space of the final states, 

= fll d 3 h d 3 k 2 \S\ 2 
T J 2 2uj (2vr) 3 2wi (27r) 3 2w 2 VT ' [ ' 

This can be expressed as 




\M\ 2 uu 1 uj2d 3 k 1 d s k26(Y^ Ui)6(Y^ K) . 



where the matrix element M is a function of energies and polarization vectors of the incoming 
and outgoing photons. 



3.1. Matrix elements 

We use the system of units in which Ti = c = m e = 1, so the fine structure constant is a = e 2 . 
The magnetic field is expressed in the units of the critical field B c . The effective Lagrangian of the 
electromagnetic field with QED corrections is given by (Berestetskii, Lifshits and Pitaevskii 1982) 

L eS = ^j o -^-xl-^cotO(vcoth V ) + l-^L\ . (22) 

Here 

and F = \ (B 1 - £ 2 ) , J = BE are the field invariants. In the low frequency approximation the 
S-matrix element can be calculated as 

"S-y— ►71+72 = (7 I / drdtVint|7i72>, 
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where the interaction operator is VS nt = L c g. 

We find that the two lowest order terms of the S-matrix correspond to a square and hexagon 
diagrams: 

S & = -* (2 " )2315 (4tt) 3 / 2 x (2vr) 4 ^ 2 5(]>>^(]>>)M 6 . (23) 

and 

S 4 = (47r) 3 / 2 i?^ 1 ^ 2( 5(^ Wi ) <5(£ fei) M 4 . (24) 

The matrix elements are 

M 6 = 48^ + 26,6 + 13C + 162? , 
M 4 = 8C + 14D, 

and denoting the by rtj, rij, e*, the index of refraction, the direction of propagation, and the 
polariztion vector of th i— th photon respectively we obtain: 

.4 = n nin 2 (n x e )*(rai x e[) z (n 2 x e 2 ) 2 

6 = n (n x eo) z e[ z e% z + ni(ni x e[) z e 0z e* 2z + n 2 (n 2 x e* 2 ) z e Qz e\ z . 

C = n (n xe ) z {[n 1 n 2 (n 1 n2)-i\(e* 1 e* 2 )-n 1 n2(n 1 e* 2 )(n 2 e* 1 )} + 
ni(ni x e*) 2 {[non 2 (non 2 ) - l](e e 2 ) - n n 2 (n e 2 )(n 2 e )} + 
n 2 (n 2 x e 2 ) 2 {[n ni(n n 2 ) - l](e e*) - n ni(n ei)(nie )} 

and 

V = e Qz (nini - n 2 n 2 )(e\ x e 2 ) + e^(n n - n 2 n 2 )(e x e 2 ) + e 2z (n n - n 1 n 1 )(e x e*) 
The matrix element can be expressed as 

S = - i J^(^) 3/2 ^i^M(27T) 4 5(J2^)5(J2k l ) (25) 

where 

M = £s Me + 4^ M <- < 26 > 

3.2. Polarization selection rules 

In the vacuum, polarization of the normal modes is linear, i.e.: 

n = (cos 9,0, sin 9) 

ey = (sin 0,0, -cos6>) (27) 
e± = (0,1,0). 
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Functions C = and T> vanish when we neglect refraction, while otherwise they are of the order of 
n — 1. Using equations ( p7| ) and neglecting refraction we find that M4 vanishes while 

M 6 (_L^_L + _L) = 48sin 3 6» 
Af 6 (±-HI + II) = 26sin 3 # 
M 6 (||->_L + ||) = 26sin 3 0, 

and Mq = for all other transitions. 

In the presence of matter the functions C and T> no longer vanish. This is due to the fact that, 
in general, the refraction coefficients differ from unity, polarization vectors are elliptical rather 
than linear, photons in the process are not exactly collinear, and there exists a small degree of 
longitudinal polarization. The polarization related effects are the strongest for high matter density 
and propagation direction along the field. 



3.3. Kinematic selection rules 

The kinematic selection rules arise from the fact that for some transitions the the energy 
conservation cannot be satisfied. The kinematic selection rules can be discussed analytically using 
the dispersion relation k = n(tu,fi)uj = (1 + /u))w. We note that the refraction coefficient 
is very close to unity, see Figures |l| and and therefore for a given [i and one can invert the 
dispersion relation to obtain 

UJ = (l-Z(k,ri + 0(e))k (28) 

We note that for the fixed initial photon the function / = ujq — oj\ — u>2 is monotonically decreasing 
with //- the cosine of the angle between the initial and one of the final photons. Thus, it suffices 
to verify that equation / = has no solution for the collinear photons, to be sure that there are no 
solutions for non-collinear photons. Using momentum conservation and equation (|28|), one obtains 
for the collinear photons: 

/collinear = -£(&(), W)) + £(h, Ho)h +£(&2,/Uo)&2 (29) 

The condition for the kinematic selection rules to be satisfied is / co iiinear > 0. 

In vacuum in the limit ^ < 1 the refraction coefficients are functions of the magnetic field 



and the dependence on photon energy to is weak, see Figure In this case condition (29) can be 
rewritten as — £o + x£,i + (1 — ^)^2 > 0, where x is a number between and 1. It is clear, that 
transitions II— >|| + _L and II— >_L + _L are not allowed. Since in the first order in — the functions 
£ increase as a function of u transitions ||— »|| + || and _L^_L + _L are also forbidden. Thus in 
vacuum only two transitions are kinematically allowed: _L^_L + || and _L— >\\ + ||. 

Combining the polarization selection rules and the kinematic selection rules we conclude, in 
agreement with Adler, that only the transition _L— >\\ + || is allowed vacuum. 
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In the presence of matter the discussion of kinematic selection rules becomes more complicated 
since the refraction indices as a function of energy are rather complicated functions. The kinematic 
selection rules become complicated when the field is strong enough so that the electron cyclotron 
resonance influences refraction significantly. Around the resonance refraction coefficients are non 
monotonic functions of energy. In this case there may be a fraction of final state space that 
becomes kinematically allowed because of the influence of plasma on refraction coefficients. 



3.4. Absorption coefficient 

In the case of propagation in vacuum we ignore dispersion and calculate the integral of 
equation ( ^Tl) 



/ 



2q,6 r q,6 ^5 

dr = ^ )3 ^| 2 x J kki(k ~~ k{) kfdkid cos 6 d(j)S(k - k\ — \k - k\\) = ^ — ^ 1-^^ 1 2 • 



Due to the selection rules discussed above the only non-vanishing photon splitting absorption 
coefficients is 

,. (1H i + n )= 4 B o si „«^(^\ 

v 11 UJ 2tt 2 30 V315/ 

in full agreement with the result obtained by Adler (1971). 

In general, when the polarization modes are not linear and we take into account the dispersion 
relation equation ( pi] ) can be written as 

2a 6 f 

IM^uJuii^kfdkidcosO'dcfi'd (oj — oj\ — 0J2) (30) 



(2vr) J 



where we wrote explicitly the variables in the system U' and uj-s are function of k and \x. Defining 
/ = ujq — u>i — u>2 , we can integrate over cos 6' and obtain 



2a 6 

(2^) 3 7o 



2tt 



dk\k\ \ dcj) bJujibj 2 \MY 



df 



(31) 



which can be evaluated numerically. The integrand in equation ([3l]) is understood to vanish 
whenever there is no solution of equation / = 0. 



4. Results 

We evaluate equation (^) numerically to obtain the photon splitting absorption coefficient for 
different polarization channels, and various photon energies, magnetic fields and plasma densities. 
In the calculation we use the refraction indices of equation (||), the polarization vectors given by 
equation (||), and calculate the matrix element using equation (^). At each integration point we 
evaluate the function / and find whether the energy conservation is satisfied. 
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We present the results in Figures H [|, [|, and ||. Each figure consists of four panels which 
show the splitting rates for four angles of propagation with respect to the magnetic field: 10°, 20°, 
30°, and 70°. Figures [| || show the splitting absorption coefficient as a function of matter density 
of a photon with the energy Tiuj = 1.5m e c 2 ; Figure || for the case B = B c and Figure || for the 
case B = 2B C . We present photon splitting absorption coefficients as functions of energy when the 
plasma density is lOOg cm -3 in for B = B c in Figure [5L and for B = B c in Figure |(| We show the 
photon splitting absorption coeficients for the energies below 2m e c 2 since above this energy the 
opacity is dominated by the single photon pair production. 

We first consider the effects of matter density on the photon splitting rates. At low densities 
the influence of matter is negligible, and we recover the vacuum case when all the polarization 
vectors are linear, and the kinematic selection rules are determined by the magnetic vacuum 
refraction. With the increase of density polarization of the normal modes becomes elliptical 
starting at the photons propagating near to the direction of the field. At small angles this effects 
are pronounced already at the density of O.lg cm -3 . However comparing Figures 4 and 5 we 
see that the dominant effect is due to the kinematic selection rules and the influence of the 
electron cyclotron resonance is crucial. This is also clearly seen in Figure 6, where a number of 
splitting channels that are forbidden in vacuum suddenly turns on for energies above the cyclotron 
resonance. 

At a low value of the magnetic field the process rapidly becomes unimportant, since the rate 
scales -B 6 . With the increasing value of the magnetic field the plasma effects start to be important 
when the electron cyclotron resonance falls right around the electron mass, i.e. the field has a 
value close to the critical field. In these case when the matter is sufficiently dense the effects of the 
electron cyclotron resonance influence the kinematic selection rules significantly thus allowing more 
photon splitting polarization channels. When the value of the magnetic field is higher than 2B C 
the electron cyclotron resonance falls above region of integration over the final states and the does 
not influence the splitting rate. Moreover with the increasing magnetic field refraction becomes 
dominated by the vacuum terms, and the effects of plasma become less and less important. 

Figures |5| and ^ show the effects of photon energy on the splitting absorption coefficient. 
When the electron cyclotron resonance effects are ignored only small frequencies and small angles 
of propagation are influenced, since it is there where the refraction coefficients and the polarization 
of the normal modes are the most influenced by plasma. However this is also the region where the 
photon splitting coefficient is the smallest, see equation (|I|). 

5. Discussion 

In this work we extend the results of Adler (1971) to the case of propagation in the magnetized 
plasma and concentrate on the case of the density typical for a neutron star atmosphere. Our 
approach is accurate for the magnetic fields up to approximately the critical field since we use 
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a low field approximation in the calculation of the matrix element. We calculate the refraction 
coefficients accurately and thus the kinematic selection rules do not suffer from this limitation. 

We have calculated the photon splitting absorption coefficient as a function of magnetic field, 
plasma density, and the photon energy and direction. We have found a region of the parameter 
space (density p > lg cm~ 3 , the magnetic field 0.1 B c < B < 2B C , propagation angles 9 < 30°) 
where the effects of plasma are the most pronounced. A part of these region is where the photon 
splitting absorption is small, so the region of importance is limited to photon energies above m e c 2 . 
We find that the the photon splitting rate is well described by the vacuum approximation in the 
remaining part of the parameter space. 

Photon splitting absorption coefficient is small when compared to other processes 
that may play a role in plasma, for example the electron scattering opacity is 
kscat = 0A(huj/m e c 2 ) 2 (B c /B) 2 (p/g cm~ 3 )cm _1 for the extraordinary mode, a value a few 
orders of magnitude higher than that for photon splitting. Thus photon splitting can play a 
significant role only in very special astrophysical cases. An example of such environment could be 
in deep layers of a neutron star atmosphere. Therefore the results presented here may apply to soft 
gamma-ray repeaters, where a large amount of energy is deposited in the crust of a neutron star. 
Photon splitting in a high density plasma may be a way of producing a large number of soft X-ray 
photons which later on escape. Our results can also be applied to the high energy radiation from 
isolated neutron stars provided that we see radiation from the surface and not the magnetopshere. 
However, the main conclusion is that surprisingly the effects due to the presence of plasma are 
important in a rather small fraction of the parameter space and the vacuum approximation can 
be used in most calculations. 
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Fig. 1. — Refraction coefficients in magnetized vacuum as a function of the magnetic field, for 
propagation angle of 90°. The upper line corresponds to the _L mode and the lower curve to the || 
mode. 
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Fig. 2. — Refraction coefficients in magnetized vacuum as a function of photon energy for three 
values of the magnetic field: B = 0.3 x B c (solid lines), B = B c (doted lines), and for B = 10 x B c 
(dashed lines). For each magnetic field the upper line corresponds to the _L mode and the lower 
curve to the || mode. The direction of propagation is perpendicular to the magnetic field. Note 
that at photon energies above m e C 2 the low energy approximation n — n(uj = 0) oc uj 2 no longer 
holds, especially for the very strong field case. 
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Fig. 3. — Photon splitting absorption coefficient as a function of density in plasma for B = B c . 
The solid line corresponds to the case 1^2 + 2, i.e. the transition allowed in vacuum, the dotted 
line represents the transition 2 — > 1 + 2, the long dashed 2 — ► 2 + 2, the short dashed dotted 
1^2 + l,and the long dashed dotted 1 — ^ 1 + 1 
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Fig. 4. — Photon splitting absorption coefficient as a function of density in plasma for B = 2 B c . 
For explanation see Figure 3. 
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Fig. 5. — Photon splitting absorption coefficient as a function of energy in plasma with density 
100 g cm~ 3 for B = B c .Fov explanation see Figure 3. 
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Fig. 6. — Photon splitting absorption coefficient as a function of energy in plasma with density 
100 g cm -3 for B = 2 B c . For explanation see Figure 3. 



